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Multiple scattering formulations have been employed for more than 30 years as a method of 
studying the quantum vacuum or Casimir interactions between distinct bodies. Here we review the 
method in the simple context of (5-function potentials, so-called semitransparent bodies. (In the 
limit of strong coupling, a semitransparent boundary becomes a Dirichlet one.) After applying the 
method to rederive the Casimir force between two semitransparent plates and the Casimir self-stress 
... on a semitransparent sphere, we obtain expressions for the Casimir energies between disjoint parallel 

OO ' semitransparent cylinders and between disjoint semitransparent spheres. Simplifications occur for 

' weak and strong coupling. In particular, after performing a power series expansion in the ratio 

, of the radii of the objects to the separation between them, we are able to sum the weak-coupling 

expansions exactly to obtain explicit closed forms for the Casimir interaction energy. The same can 
' be done for the interaction of a weak-coupling sphere or cylinder with a Dirichlet plane. We show 

^ , that the proximity force approximation (PFA), which becomes the proximity force theorem when 

■ the objects are almost touching, is very poor for finite separations. 

^ ' PACS numbers: 03.70.-|-k, 03.65.Nk, ll.80.Et, ll.80.La 



p . I. INTRODUCTION 

I) : 

Recently, there has been a flurry of papers concerning "exact" methods of calculating Casimir energies or forces 
[ between arbitrary distinct bodies. Most notable is the recent paper by Emig, Graham, JafFe, and Kardar [l|. (Details, 
^SJ . applied to a scalar field, are sup^hed in Ref. 0. See also Refs. 0, Q-) Precursors include an earlypaper of Renne [1], 
^ ■ rederiving the Lifshitz formula [6] in this way, the famous papers of Balian and Duplantier 0, S , work of Kenneth 
, and Klich |1Q| ba sed on the Lippmann-Schwinger formulation of scattering theory [llj , the papers by Bulgac, Marierski, 
and Wirzba ITsl [T3 | , who use the modified Krein formula [iBl , and by Bordag [l^ fT7| , who derives his results 
from a path integral formulation. Dalvit et al. p^.[l9j use the argument principle to calculate the interaction between 
conducting cylinders with parallel axes. See also Reynaud et al. '20] and references therein. 

In fact, Emig and earlier collaborators [2l|, [l^, have published a series of papers, using closely related methods 
to calculate numerically forces between distinct bodies, starting from periodically deformed ones. Strong deviation 
from the proximity force approximation (PFA) is seen for when the distance between the bodies is large compared 
to their radii of curvature. Bordag has precisely quantified the first correction to the PFA both for a cylinder 
and a sphere near a plane. As Gies and Klingmiiller note [l^], 1% deviations from the PFA occur when the ratio 
k> ] of the distance between the cylinder and the plate to the radius of the cylinder exceeds 0.01. We will not discuss 
5_j ' the worldline method of Gies and collaborators [25', ^2^, |23| further, as that method lies rather outside our discussion 
^ here. Similar remarks apply to the work of Capasso et al. [1^, who calculate forces from stress tensors using the 
familiar construction of the stress tensor in terms of Green's dyadics [1^ Is^l ■ They use a numerical finite-difference 
engineering method. 

It is clear, then, with the exception of these last two methods, these approaches are fundamentally equivalent. We 
will refer to all of the former methods as multiple scattering techniques. We will now proceed to state the formulation 
in a simple, straightforward way, and apply it to various situations, all characterized by 5-function potentials. (A 
preliminary version of some of our results has already appeared (sH-) 
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II. FORMALISM 



We begin by noting that the multiple-scattering formalism may be derived from the general formula for Casimir 
energies (for simplicity here we restrict attention to a massless scalar field) [32| 

S-^TrlnG, (2.1) 

where r is the "infinite" time that the configuration exists, and G is the Green's function in the presence of a potential 
V satisfying (matrix notation) 

[-d^ + V)G=l, (2.2) 

subject to some boundary conditions at infinity. (For example, we can use causal or Feynman boundary conditions, or, 
alternatively, retarded Green's functions.) In Appendix [Al we give a heuristic derivation of this fundamental formula. 

The above formula for the Casimir energy is defined up to an infinite constant, which can be at least partially 
compensated by inserting a factor as do Kenneth and Klich [lol |: 

S=^TrlnGGo\ (2.3) 

Here Go satisfies, with the same boundary conditions as G, the free equation 

- a^Go = 1. (2.4) 

Now we define the T- matrix (note that our definition of T differs by a factor of 2 from that in Ref. [l|) 

T ^ S -1 = V{1 + GqV)-^. (2.5) 

We then follow standard scattering theory as reviewed in Kenneth and Klich [T3|. (Note that there seem to be 
some sign and ordering errors in that reference.) The Green's function can be alternatively written as 

G = Go - GoTGo = ^ , \ Go = ^-^TGo, (2.6) 
i + Cto y 

which results in two formulas for the Casimir energy 

E = — Trln (2.7a) 

2r 1 + GoV ^ ' 



— TrlnV^-^T. (2.7b) 



If the potential has two disjoint parts, 
it is easy to show that 



V = Vi + V2, (2.8) 
T = (l/i + V2){1 - GoTi)(l - GoTiGor2)-i(l - GqTs), (2.9) 



where 

T,^V,{\ + GoV,)-\ z = l,2. (2.10) 

Thus, we can write the general expression for the interaction between the two bodies (potentials) in two alternative 
forms: 

Ei2 = -:^Trln(l - GoTiGoTa) (2.11a) 
It 

= __LTrln(l-FiGiV2G2), (2.11b) 

where 

G, = (l + GoV^O"'Go, z=l,2. (2.12) 

The first form is exactly that given by Emig et al. [l|, and by Kenneth and Klich JXjj, while the latter is actually 
easily used if we know the individual Green's functions. (The effort involved in calculating with either is identical.) 
In fact, the general form (|2.11ap was recognized earlier and applied to planar geometries by Maia Neto, Lambrecht, 
and Reynaud [l^, [H, [s^l ■ In fact, Renne [9| essentially used Eq. (j2.11bp to derive the Lifshitz formula in 1971. 
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III. CASIMIR INTERACTION BETWEEN 5-PLATES 

We now use the second formula above (|2.11bp to calculate the Casimir energy between two parallel semitransparent 
plates, with potential 

V = Xi6{z- zi) + \25{z- Z2), (3.1) 

where the dimension of Ai is L^^. The free reduced Green's function is (where we have performed the evident Fourier 
transforms in time and the transverse directions) 

g^{z,z')^^e--\'-''\, K^ = Q^ + e. (3.2) 

Here k = is the transverse momentum, and C, — —iw is the Euclidean frequency. The Green's function associated 
with a single (5-function potential is 



Then the energy/area is 



S = I dC I d^k I dz\n{l~ A){z,z), (3.4) 



where, in virtue of the (5-function potentials {a — \z2 — zi\) 

-'■"■"■-"■'^'-^''<5(z-zi). (3.5) 



Ai + 2k A2 + 2k 
We expand the logarithm according to 



00 ^ 



ln(l-A) = -^— . (3.6) 



.=1 ' 

For example, the leading term is easily seen to be 

^(2) _ [ dC <Pk _2Ka _ ^1^2 I . 

~ 16^3 7 ^ " 32^2^' ^'^■'^ 

which uses the change to polar coordinates, 

dCd^k = dKK^ dVL. (3.8) 
In general, it is easy to check that, because A{z, z') factorizes here, A{z, z') — B{z)C{z'), Tr^" — (Tr^)", or 

Trln(l-A) = ln(l-TrA), (3.9) 
so the Casimir interaction between the two semitransparent plates is 



dKK^ln 1-- ^e-"°- —e-'"'], (3.10) 

47r2 L V Ai + 2k A2 + 2k ' ' ^ ' 



which is exactly the well-known result [3£ 



IV. CASIMIR SELF-ENERGY FOR A SINGLE SEMITRANSPARENT SPHERE 



Before we embark on new calculations, let us also confirm the known result for the self-stress on a single sphere 
of radius a using this formalism. (This demonstrates, as did the rederivation of the Boyer result [3(| by Balian and 
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Duplantier 's'l, that the multiple scattering method is equally applicable to the calculation of self-energies.) We start 
from the general formula (|2.7ap . where 

F(r,r') = A(5(r-a)5(r-r'). (4.1) 
We use the Fourier representation for the propagator in Euclidean space, 

e-|Clk-r'| r ffj^ gik-(r-r') 

as well as the partial wave expansion of the plane wave 

e^kr ^J24Wjiikr)YUr)Yi:ri^)- (4-3) 

Im 

Then, from the orthonormahty of the spherical harmonics, 

(4.4) 



(4.5) 



we obtain the representation 



2 . dh 

Go(r,r') - -E/, ^j:j:^J/(A;r)jHfc^')n„.(f)l^™(f')- 



Now we combine the representation for the free Green's function with the spherical potential (|4.ip to obtain 

When this, or powers of this, is traced (that is, r and r' are set equal, and integrated over), we obtain a poorly defined 
expression; to regulate this, we assume r ^ a, for example, r < a. (This is a type of point-split regulation.) Then, 
because 

jiika) = i (/»«(fca) + h['\ka)) = ^ [hl'^ka) + (-l)'/i[^)(-fca)) , (4-7) 
while ji{kr) = (— l)'j/(— /cr), we see that the k integration in Eq. (|4.6p can be evaluated as"'^ 

^j^Ji(fc«)jKfc'') = -A'/+i/2(lCl«)/i+i/2(ICk), r<a. (4.8) 
Thus, it is easily seen that an arbitrary power of GqV has trace 

Tr {G,Vr = {Xar ^ (^/+i/2(IC|a)/;+i/2(IC|a))" , (4.9) 

Im 

and that therefore the total self-energy of the semitransparent sphere is given by the well-known expression [s^, [s^ 

E^—y^{2l + l) dx\n{l + XaIi+^/2{x)Ki+y2ix)) , x ^ \C\a. (4.10) 
^7ra Jo 

Actually, a slightly different form involving integration by parts was given in Refs. [s^ |40| , which results in the energy 
being finite though order . In order A"^ there is a divergence which is associated with surface energy [41i| . 



Of course, this result is the immediate consequence of the usual partial wave expansion 

Go (r, r') = ikJ2 h {kr<)h^^^ {kr> )Yi^ {r)Y;^ (f k = \uj\. 

Im 

The point of our slightly more elaborate approach hero is that it generalizes to the corresponding two-body case — see Eq. | |5.7t . 
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V. 2 + 1 SPATIAL GEOMETRIES 



We now proceed to apply this method to the interaction between bodies, which leads, for example, as Emig et al. 
[IliB point out, to a multipole expansion. In this section we illustrate this idea with a 2 + 1 dimensional version, which 
allows us to describe, for example, cylinders with parallel axes. We seek an expansion of the free Green's function for 
R = R_L entirely in the x-y plane. 



Go(R + r'-r) = 



|r-R-r 



47rr-R-r' 



dk 



2tt 



where the reduced Green's function is 

go{ri_ - Ri - t'l) = 



(d^fci) e 



~ik±-R± ik.±-{r±-r'^) 



(27r)2 fci + fc2 + C 



(5.1) 



(5.2) 



As long as the two potentials do not overlap, so that we have rj^ — R^ — 7^ 0, we can write an expansion in terms 
of modified Bessel functions: 
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(r^ - Rj_ - rl) = ^ /™(^r)e""^/™,(Acr')e-™'^'ff^,„,(^i?), = fc^ ^ ^2 



(5.3) 



By Fourier transforming, and using the definition of the Bessel function 



^—iriKp ^ikr cos <; 



we easily find 



271 Jq k? + k 



2tt 



dkk i^^-^^ Jm{kr)Jm'{kr') 



(5.4) 



(5.5) 



which is in fact independent of r, r' . 

As in the previous section, the k integral here can actually be evaluated as a contour integral, as Bordag noted 
No point-splitting is required here, because the bodies are non-overlapping, so r/R, r' / R < 1. We write the dominant 
Bessel function in terms of Hankel functions, 



J n 



(1) 



and then we can carry out the integral over k by closing the contour in the upper half plane. We are left with 
dx X 



9 , 9 ^m—m' 





ix)J„^{xr/R)J,n'{xr'/R) = (-1)™ K^_,^,{y)I^{yr / R)I„,,{yr' / R), 



and therefore the reduced Green's function has the simple form 

f— IV"' 

ZTT 



(5.6) 



(5.7) 



(5.8) 



Thus we can derive an expression for the interaction energy per unit length between two bodies, in terms of discrete 
matrices. 



e^^ = ^JdCdk, Indet (1 ~ ~g%~9'^h) , 
where T denotes transpose, and where the T matrix elements are given by 

drrd(j) [ dr' r' d(j)'Im{Kr)e-'""f'I,n,{Kr')e'"'''l''T{r,(j);r\(l)'). 



(5.9) 



(5.10) 
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A. Interaction between semitransparent cylinders 

Consider, as an example, two parallel semitransparent cylinders, of radii a and 6, respectively, lying outside each 
other, described by the potentials 

Vi^\i5{r-a), V2 = \25{r' -h), (5.11) 

with the separation between the centers R satisfying R > a + b. It is easy to work out the scattering matrix in this 
situation, 

Ti = Vi- ViGoVi + ViGoViGoVi (5.12) 
so the matrix element is easily seen to be 

(il)mm' — 2lTXiaSmm' . , . T^T \ I \- (5.13) 

Again, we used here the regularized integral^ 

dkk 







Thus the Casimir energy per unit length is 



Jm{kr)J,n{kr') = K^{nr')I^{Kr), r < r' . (5.14) 



1 

€=— y dKKtrln(l - A), (5.15) 



where 



in terms of the matrices 



A = B{a)B{b), (5.16) 
B„,^,{a) = XiaPA^a) ^^^^^^ 

1 + Aialm' [Ka)K,n'{K.a) 

B. Interaction between cylinder and plane 

As a check, let us rederive the result derived by Bordag for a cylinder in front of a Dirichlet plane perpendicular 
to the X axis. We start from the interaction (|2.11ap written in terms of G2, the deviation from the free Green's 
function induced by a single potential, 

G2 ~ G2 — Go = —GoT2Go, (5.18) 

so the interaction energy has the form 

£; = -;^Trln(l + riG2). (5.19) 
It 

When the second body is a Dirichlet plane, G may be found by the method of images, with the origin taken at the 
center of the cylinder, 

G(r,r') = -Go(r,r'), f' ^ {R- x\y\z'), (5.20) 



^ Again, this is equivalent to the use of the two-dimensional Green's function 

00 

Ho(kP)= i'"e-''"^'j„{fcp')e™'^J^m'(M, P'<P> 

m — — 00 



where P = -sj ^ p''^ — 2pp' cos{</) - (p'). 
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where R is the distance between the center of the cyhnder and its image at Rj^, that is, R/2 is the distance between 
the center of the cyhnder and the plane. (We keep R here, rather than R/2 — D, because of the close connection to 
the two cylinder case.) Now we encounter the 2-dimensional Green's function 

giri_ + rY - R^) - ^ IraiKr)IrrAi^r')e'"'^e'"''^' gmm'{KR), (5.21) 

mm' 

(because the cylinder has y —y reflection symmetry) where the argument given above yields 

gmm'{liR) = ^Km+m'{nR)- (5.22) 

Thus the interaction between the semitransparent cylinder and a Dirichlet plane is 

1 r°° 

— / KdKtr ln(l - B(a)), (5.23) 
47r J a 

where B{a) is given by Eq. (|5.17p . In the strong-coupling limit this result agrees with that given by Bordag, because 

tr = tr B^ = ] . K^+m' {nR)Im' [no). (5.24) 

C. Weak-coupling 

In weak coupling, the formula (j5.15p for the interaction energy between two cylinders is 

\ X ab °° f°° 

(S^-^^ E J dxxKl^„A^)ll{^alR)llA^b/R). (5.25) 
Similarly, the energy of interaction between a weakly-coupled cylinder and a Dirichlet plane is from Eq. (|5.23p 

6: = -^ E j dxx K2m{x)ll{xa/R). (5.26) 



m— — oo 



D. Power series expansion 

It is straightforward to develop a power series expansion for the interaction between weakly-coupled semitranspar- 
ent cylinders. One merely exploits the small argument expansion for the modified Bessel functions Im{xa/R) and 
I^.{xb/R): 



(rp \ 2\m\ , / ■T \ 2ri 

I) E^l™l."(f) ' (5.27) 

n=0 



where the coefficients Zm,n are 



E 

fc=0 



fc! {n - fc)! r(fc + m + 1) T{n ~ k + m + 1) 



22(m+n) r(m + n+ i) 
~ y/n nl (2m -|- n)\ T{m + n+l)' 

The Casimir energy per unit length (|5.25p is now given as 

oc oo oo oo I / 7 \ 2|m'| / , \ 2n' 



(5.28) 



_ AiaA26 
47ri?2 



y E E EES) ^H.n(l^) (^) 

^ m—~oo m' — — oo n— n'— 



(5.29) 
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Reordering terms gives a more compact formula 



iaX2b ^ ^ /a\2(|m|+n) 

m— — oo 7n' — — oo n— n'— 



7i ,1 



^ X 2(|m'|+n') 



|m|4-|?Ti'|+'n'+^T''-'m,+'m,' ' 



where the two index symbol Jp^q represents the integral over x, which evaluates to 



7 — 2 



^p\Tip + q+l)T{p-q+l) 
22p+2r(p+|) 



(5.30) 



(5.31) 



In order to simplify the power series expansion in terms of and ^ we need to reorder the m sums so that only 
non-negative values of m appear. There are several ways to break up the m sums; one of them is to decompose the 
sum into the m = m' = term, the m, m' same sign terms, and the m, m' different sign terms, giving 



AiaA26 



OC OC 



2n 



n=0 n'=0 

OO OO OO OO 



2n' 



J, 



n-\-n' .0 



+2EEEE^-Gy 

m— 1 m'— n— n'— 

OO OO OO OO 

+2EEEE^-(i) 



a \ 2(m+n) 

/ a \ 2(m+n) 



2(m'+n') 



m— m' — 1 n— n'— 



ra-\-ra' +n+n' ,m+m' 
2{m'+n') 



J, 



(5.32) 



It is now possible to combine the multiple infinite power series into a single infinite power series, where each term is 
given by (possible multiple) finite sum(s). In this case we get an amazingly simple result 



\1a\2b 1 / a 



n=0 



2n 



where fi = b/a, and where by inspection we identify the binomial coefhcients 

I'M -til 



k=0 



,2k 



(5.33) 



(5.34) 



Remarkably, it is possible to perform the sums [43 | , so we obtain the following closed form for the interaction between 
two weakly-coupled cylinders: 



AiaA26 



1 - 



a + b 
R 



1 - 



~r' 



-1/2 



(5.35) 



We note that in the limit R — a — b = d ^ 0, d being the distance between the closest points on the two cylinders, we 
recover the proximity force theorem in this case (jB4|, 



Uid) 



A1A2 2ab 1 



327r V R d'^l'- 



■, d<s: a, b. 



(5.36) 



In Figs.[TH2]we compare the exact energy (I5.35|) with the proximity force approximation (|5.36p . Evidently, the former 
approaches the latter when the sum of the radii a -I- 6 of the cylinders approaches the distance R between their centers. 
The rate of approach is linear (with slope 3/2) for the equal radius case, but with slope b^ /4a^ when a <^ b. More 
precisely, the ratio of the exact energy to the PFA is 



[/ ~ Afi R 



1 - 



R'^-aR + a^ d 
Aa{R~a) R' 



(5.37) 



This correction to the PFA is derived by another method in Appendix [Cl The reader should note that the the PFA 
is actually only defined in the limit d — s- 0, so the functional form away from that point is ambiguous. Corrections to 
the PFA depend upon the specific form assumed for U(d). 
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FIG. 1: Plotted is the ratio of the exact interaction energy (|5.35|) of two weakly-coupled cylinders to the proximity force 
approximation (|5.36|l as a function of the cylinder radius a for a = b. 




0.000 0.002 0.004 0.006 0.008 0.010 
a/R 

FIG. 2: Plotted is the ratio of the exact interaction energy l|5.35[) of two weakly-coupled cylinders to the proximity force 
approximation (|5.36|l as a function of the cylinder radius a for b/a — 99. 



E. Exact result for interaction between plane and cylinder 



In exactly the same way, starting from Eq. (j5.26p . we can obtain a closed-form result for the interaction energy 
between a Dirichlet plane and a weakly-coupled cylinder of radius a separated by a distance R/2. The result is again 
quite simple: 



Xa 
'47ri?2 



1 - 



-3/2 



In the limit as d ^ 0, this agrees with the PFA: 



Uid) 



X V2a 

'64^^372' 



(5.38) 



(5.39) 



Note again that this form is ambiguous: the proximity force theorem is equally well satisfied if we replace a by -R/2, 
for example, in U{d). The comparison between this PFA and the exact result (|5.38|) is given in Fig. [31 



F. Strong coupling (Dirichlet) limit 

The interaction between Dirichlet cylinders is given by Eq. (I5.15[) in the limit Ai = A2 00, that is 

1 r°° 

(B= 7T dxxtrln(l - A), 



(5.40a) 
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1.0 

0.8 ; 

;d 0.6 - 




0.0 0.1 0.2 0.3 0.4 0.5 



a/R 

FIG. 3: Plotted is the ratio of the exact interaction energy (|5.38|) of a weakly-coupled cylinder above a Dirichlet plane to the 
proximity force approximation (|5.39p as a function of the cylinder radius a. 



^ Km" [xa/R) Km' (xb/R) 



where 



Here the trace of the logarithm can be interpreted as in Eq. (|3.6p . 

Because it no longer appears possible to obtain a closed-form solution, we want to verify analytically that as the 
surfaces of the two cylinders nearly touch, we recover the result of the proximity force theorem. We use a variation of 
the scheme explained by Bordag for a cylinder next to a plane (l6j . [The analysis is a bit simpler in the weak-coupling 
case, which leads to Eq. (|5.36|) . See Appendix [Cl] First we replace the products of Bessel functions in A by their 
leading uniform asymptotic approximants for all m's large: 

/ / X 2\ -1/4 / , , 2\ -1/4 

1 1 1 fx //-X 



Bmm"ia)Bm"rn'{b) ^ _= — == M + „ 1 + , (5.41) 



where the exponent is 

X = (m + to")?/ ( ) + (^"' + ™")'? ( "TT — 7j) ~ 2"^"^' ("TTh) ~ '^'^''^ ( ~T5 

\m + ni" J \m' + m" J \m"R/ \m'R 

in terms of 



1=1 



This normalization is chosen so at the critical point where x — ^ for a + b = R^ 



(5.42) 



riiz) = t-'+ln-^, V'{z) = -^, ry"(z) = -4, (5.43) 
1 + < ^ zt 

and 

t = (l + z2)-i/2. (5.44) 
We write the trace of the sth power of A as (summed on repeated indices) 

{A )mimi = -Bmim'j (o)-Bm'^m2 (^)-^m2m2 ('^)^™2™3 ' ' ' Bm^m'^^O^Brn'^miib) . (5.45) 

We rescale variables in terms of a large variable M and relatively small variables: 

m- = Mar, m = MP,, (5.46) 
where without loss of generality we take only 2s — 1 of the a's and /3's as independent: 



^(«^+A) = s. (5.47) 



a. = | A = l-|' Vz. (5.48) 
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Away from this point, we consider fluctuations, 



with the constraint 



Tl It 



= 0. 



(5.49) 



(5.50) 



The Jacobian of this transformation is sM^" ^ . 

Now, we expand the exponent in tr^**, to first order in d = R — a — b, and to second order in cti, [3i. The result is 



2Msd fR , \ 



tR 



^ (A + A+i) 

z It — a 



Mt a 
~ R-a 



(5.51) 



The ai terms lead to trivial Gaussian integrals. The difficulty with the quadratic Pi terms is that only s — 1 of the 
differences are independent. But, in view of the constraint (j5.50p there are only s — 1 independent (3i variables. In 
fact, it is easy to check that 



s-l 



E(A-A+i)^ = E 



i + l 



1 



(5.52) 



which now enables us to perform each successive (3i — integration. The Jacobian of the transformation to the 
difference variables Ui — f3i ~ /3i+i, « = 1, . . . , s — 1, is 1/s. Thus, we can immediately write down 



-1 poo ^ ,s foo j\,r2s+l 

AnR^ Jo ^ Wo i^^MY 



^-2Msd/tR 



-I s s-l 



i±i</4(_R-a) 



47ri?2 



dz z ■ 



dM M 



25 + 1 



-2Msd/tR 



{2nMy 



■na 


s/2 


'ATi(R-a)' 




{R~-a)Mt 




Mta 





(s-l)/2 



'1/2 



^2a{R-a)Tr^ / R 
3840i?3 I 7 



5/2 



(5.53) 



which is exactly the result expected from the proximity force theorem, according to Eq. (|B5p . 

We will forego further discussion of strong coupling, and presentation of numerical results, for these have been 
extensively discussed in several recent papers, especially in Ref. 



VI. 3-DIMENSIONAL FORMALISM 



The three-dimensional formalism is very similar. In this case, the free Green's function has the representation 

Go(R + r'-r)= E J/(«ICk)jH*|Ck')i"/:„(f)W(f')5w™'(R). (6.1) 



The reduced Green's function can be written in the form 



5L,w(R) = (47r)2* 



(dk) 



ik-R 



ji{kr)ji'{kr') 



(2^)3 k^ + e ]imr)]i'mr') 



-ii™(k)r;„,(k). 



(6.2) 



Now we use the plane-wave expansion ()4.3p once again, this time for e'"^'^, so now we encounter something new, an 
integral over three spherical harmonics. 



dkrz,„(k)y;„,,(k)y;,,„„(k) = c, 



Im ,1' m' ,1" m" ; 



(6.3) 
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where 

Cira,vm',i"m" = (-1) y ^ Q)\mm' m")- ^^'^^ 

The three-j symbols (Wigner coefficients) here vanish unless l + V + I" is even. This fact is crucial, since because of it 
we can follow the previous method of writing ji" (kR) in terms of Hankel functions of the first and second kind, using 
the reflection property of the latter, 

h^jf^kR) = {-lfh^^){-kR), (6.5) 

and then extending the k integral over the entire real axis to a contour integral closed in the upper half plane. The 
residue theorem then supplies the result for the reduced Green's function'^ 

= ^^^''''\[% E Ci^.v^'r^n"Ki„+„^{\C\R)Yi„^.,[A). (6.6) 

I" in" 

A. Casimir interaction between semitransparent spheres 

For the case of two semitransparent spheres that are totally outside each other, 

Vi{r)^\iS{r~a), V2(r') = A2(5(r' - 6), (6.7) 
in terms of spherical coordinates centered on each sphere, it is again very easy to calculate the scattering matrices, 

r,(r, r') = ^5(r - a)5{r' - a)^^ J'^^'i^rif nn V (^.8) 

a ^ 1 + MaKi+if2[\C\a)Ii+i/2{\C\a) 

and then the harmonic transform is very similar to that seen in Eq. (|5.13p . (fc = 

{h)im,i'm' = / idr)idr')ji{kr)Yil^ir)ji,{kr')Yi,„,,{r')T,{r,r') 



-A A ^f+i/2(K|g) 

" " ""'^ 1) 2|C| l + Aia/f,+i/2(|C|a)W(IC|a)- ^'''^ 

Let us suppose that the two spheres lie along the z-axis, that is, R = Rz. Then we can simplify the expression for 
the energy somewhat by using Yim{0 = 0) = Smox/C^l + l)/47r. The formula for the energy of interaction becomes 

^=^1 dCtYln{l-A), (6.10) 

Alm,l'm' = Sm,m'''^Bii»rn,{a)Bl"l'm,{b) (6-11) 



where the matrix 



I" 



is given in terms of the quantities 



^"""^") = 7fS^ V(2/ + l)(2/' + l)E(2^ +l)(0 jU--0j l + W.,,MC«)W(Ca) - 

(6.12) 

Note that the phase always cancels in the trace in Eq. (|6.10[) . For strong coupling, this result reduces to that found 
by Bulgac, Wirzba et al. [1^. [l4l for Dirichlet spheres, and recently generalized by Emig et al. ^] for Robin boundary 
conditions. (See also Ref. [43l|.) 



This differs by a (conventional) factor of |c^| from the quantity Ui^ii^i defined by Emig et al. Q|. 
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B. Weak coupling 

For weak coupling, a major simplification results because of the orthogonality property, 

T._^ ( m -m ) ( m -m ^0 ) = 2rTT' ^ " ^' ' ^^'^^^ 
Then the formula for the energy of interaction between the two spheres is 



E ^ - 



\ \ hi pOO 7 / 1 If Iff \ ^ 

'-htj, ^|^/^^ + ^^^^^' + ^^^^^" + ^Ko j Kl.^„^{x)ll,/,{xalR)ll^,/,{xb/R). (6.14) 



There is no infrared divergence because for small x the product of Bessel functions goes likex2('+''-'")+i, and?" < l+V . 

As with the cylinders, we expand the modified Bessel functions of the first kind in power series in a/R,b/R < 1. 
This expansion yields the infinite series 



Ai aMo ab ' ' ' 

E 



v-^ 1 v-^ /a\2(n-m)/6\ 

4.^ ^^E^TTT^^-fe) [r] ' ^'-''^ 

n=0 m=0 ^ 

where by inspection of the first several Dn.m coefficients we can identify them as 

= I l) ' (^-l^) 
and now we can immediately sum the expression (|6.15p for the Casimir interaction energy to give the closed form 

^=4-|r In f A 2 ■ (6-17) 



167ri? 



Again, when d = R — a — b <^ a,b, the proximity force theorem (jB9p is reproduced: 

t/(d) ^ MH^ln(d/i?), d<a,6. (6.18) 
lonR 

However, as Figs. |H [S] demonstrate, the approach is not very smooth, even for equal-sized spheres. The ratio of the 
energy to the PFA is 

E ^ , ln[(l + /i)V2M] 

U='+ \nd/R ' (^-^^^ 

for 6/a — /i. Truncating the power series (j6.15p at n = 100 would only begin to show the approach to the proximity 
force theorem limit. The error in using the PFA between spheres can be very substantial. 

Again we will forego discussion of the strong-coupling (Dirichlet) limit here because of the extensive discussion 
already in the literature [1, [H, [l3| . 



C. Exact result for interaction between plane and sphere 

In just the way indicated above, we can obtain a closed-form result for the interaction energy between a weakly- 
coupled sphere and a Dirichlet plane. Using the simplification that 

we can write the interaction energy in the form 

Xa 



E : 



2ttR 



/ dxj],/ — (2; + l)Xi/2(a;)/f+i/2(x(a/i?)). (6.21) 
Jo V 
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1.0 F 




FIG. 4: Plotted is the ratio of the exact interaction energy (|6.17|l of two weakly-coupled spheres to the proximity force 
approximation (|6.18p as a function of the sphere radius a for a = b. Shown also by a dashed line is the power series expansion 
(|6.15p . truncated at n = 100, indicating that it is necessary to include very high powers to capture the proximity force limit. 



1.0 P 
0.8 - 




0.000 0.005 0.010 0.015 0.020 
a/R 

FIG. 5: Plotted is the ratio of the exact interaction energy (|6.17|l of two weakly-coupled spheres to the proximity force 
approximation (|6.18|l as a function of the sphere radius a for b/a — 49. 



Then in terms of R/2 as the distance between the center of the sphere and the plane, the exact interaction energy is 

E^-^flf^L—-, (6.22) 
2tt\rJ l-(2a/i?)2' ^ ^ 

which as a ^ R/2 reproduces the proximity force Umit, contained in the (ambiguously defined) PFA formula 

C/ = -A«. (6.23) 
Stt a 

The exact energy and this PFA approximation are compared in Fig. \6\ 



VII. COMMENTS AND PROGNOSIS 

Although the multiple scattering methods date back to the 1970s ^ some remarkable new results have been 
obtained. Here we have given a perhaps simpler and more transparent derivation of the procedure than in Refs. [U, 0]- 
For example, because we have approached the problem from a general field theoretic viewpoint, we see that the 
"translation matrix" introduced there is nothing other than the free Green's function. Our approach yields the 
general form first, and the multipole expansion as a derived consequence, not the other way around. We apply this 
multiple scattering method to obtain new results for the interaction between semitransparent cylinders and spheres, 
and we have analytically demonstrated the approach to the proximity force theorem. Most remarkably, we have derived 
explicit, very simple, closed-form expressions for the interaction between weakly coupled cylinders and between weakly 
coupled spheres, as well as between weakly-coupled cylinders or spheres and Dirichlet planes. These explicit results 
demonstrate the profound limitation of the proximity force approximation, which has been under serious attack for 
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5_ 




FIG. 6: Plotted is the ratio of the exact interaction energy (|6.22|l of a weakly-coupled sphere above a Dirichlet plane to the 
proximity force approximation (|6.23p as a function of the sphere radius a. 



some time [4J, |45[. We hope that these developments will lead to improved conceptual understanding, and to better 
comparison with experiment, when they are extended to realistic materials. 
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APPENDIX A: DERIVATION OF VACUUM ENERGY FORMULA 



Following Schwinger [32l | we start from the vacuum amplitude in terms of sources, 

(0+|0_)^^' = e'^I^l, W[K] = i / {dx){dx')K{x)G{x,x')K{x'). 



(Al) 



Here G is the Green's function in the presence of some background potential. From this the effective field is 

(j){x) ^ I {dx')G{x,x')K{x'). (A2) 



If the geometry of the region is altered slightly, as through moving one of the bounding surfaces, the vacuum amplitude 
is altered: 



5W[K] ^^J {dx)idx')K{x)6G{x,x')K{x') ^ -\ j {dx)idx')(l){x)6G-\x,x')(l){x') 
which uses the fact that 

GG-^ = 1. 

Upon comparison of Eq. (jA3p with the two particle emission term in 

iW[K] ^ I J(dx)K{x)4,{x)+i J(dx)C ^ I 



i / {dx)K{x)(j){x) 



(A3) 



(A4) 



(A5) 



we deduce from the coefficient of (j){x)(j){x') that the effective two-particle source due to a geometry modification is 



iK{x)K{x') 



= -6G-\x,x'). 



(A6) 



eff 
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Thus the change in the generating functional is obtained by inserting Eq. (|A6p into Eq. (|Aip . 

SW ^^J {dx){dx')G{x,x')5G-^{x\x) = J {dx){dx')5G{x,x')G-^{x' ,x), (A7) 

which gives the change in the action under an alteration of the Green's function. From this, in matrix notation 

5W = ~-l-5TY\nG^ E=^TY\nG, (A8) 
2 2t 

because for a static configuration W — —Et, which is our starting point, Eq. (|2.ip . 

There are of course many other derivations of this famous result. For example, one can derive it rather simply 
on the basis of Schwinger's quantum action principle. It may also be worth noting that it is formally equivalent to 
another familiar representation for the quantum vacuum energy 

E = -^ r pu;'Trg, (A9) 

(for example, see Ref. [1^). Here, the Fourier transform of the Green's function appears, 

Gix,x') ^ r (*-*')t;(r,r';a;). (AlO) 

In terms of Q, Eq. (|A8p can be written in the form 

/OO J 
f^Trlng. (All) 

Because 

G'^g = 1, G'^ = -u? - + 1/, (A12) 
we see that when Eq. (|A11[) is integrated by parts, and surface terms are ignored, we immediately recover Eq. (|A9[) . 



APPENDIX B: PROXIMITY FORCE APPROXIMATION 



In this Appendix we derive the proximity force approximation (PFA) for the energy of interaction between two 
semitransparent cylinders, or two semitransparent spheres, either in the strong or weak coupling regimes. This 
approximation, relating the force between nonplanar surfaces in terms of the forces between parallel plane surfaces, 
was first introduced in 1934 by Derjaguin [4^ . This approximation is only valid when the separation between the 
bodies is very small compared to their sizes. It is now well established that the approximation cannot be extended 
beyond that limit, and that 1% errors occur if the PFA is applied when the ratio of the separation to the radius of 
curvature of the bodies is of order 1%. Fortunately, current experiments have not exceeded this limit. That should 
change in the near future, which is one reason the new numerical calculations are of importance. In fact, we have 
found that in general the errors in using the PFA may be much larger than indicated above. We concur with Bordag 
that while the proximity force theorem is exact at zero separation, any approximation based on extrapolation away 
from that point is subject to uncontrollable errors. 

Consider first two parallel cylinders, of radius a and 6, with their centers separated by a distance i? > a + 5. The 
distance of closest approach of the cylinders is d — R — a — b. The PFA consists of assuming that the energy between 
the two bodies is the sum of the energies between small parallel plane elements at the same height along the surfaces, 
that is, in polar coordinates at 9 relative to the center of cylinder a and at 9' relative to the center of cylinder b, 
where as seen in Fig. [71 

a sin 6* = 6 sin 61'. (Bl) 

Because d is much smaller than either a or b, only small values of 9 actually contribute, and the energy of interaction 
U{d) between the surface may be expressed in terms of the energy per area £{h) for the corresponding parallel plate 
problem, with separation distance h: 



U{d) = / a(i6'£:[d + a(l-cos6') + 5(l-cos6'')]- 



(B2) 
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FIG. 7: Geometry of two cylinders (or two spheres) with radii a and b, respectively, and distances between their centers of 
R > a + b. The proximity force approximation applies when the distance of closest approach d — R — a — b <^ a,b. The 
approximation consists in assuming that the interaction is dominated by the interaction of adjacent surface elements, as shown. 



Here, for weak coupling [see Eq. (|3.7p ]. 

Because is small, the PFA energy per length is 

A1A2 a 



Ai A2 



U{d) = - 



327r2 d 



de 



1 



A1A2 2ab 1 



327r V R di/2- 



(B3) 



(B4) 



To obtain the corresponding result for strong coupling, we merely replace £{h) = —77^/(1440/1"^), and a similar 
calculation yields 



U{d) 



l2ab 1 



3840 V R d^/^' 



d a.b. 



(B5) 



It is easy to reproduce the result given by Bordag [16] for a cylinder in front of a plane. For the strong coupling 
(Dirichlet) case we simply take the result (jB5[) and regard b as much larger than a, and obtain 



U{d) = - 



1920\/2d5/2^ 



d<^a. 



(B6) 



For a weakly coupled cylinder in front of a Dirichlet plane, we start from the corresponding interaction between two 
such planes, £{h) — — A/(327r^/i^), which leads to 



Uid) 



A (2a)V2 
647r d3/2 



(B7) 



For nearly touching spheres the calculation goes just the same way. The result, for strong coupling (Dirichlet 
boundary conditions), for the PFA energy is 



U{d) = - 



ab 1 



d <^ a,b, 



1440 R d2 

while in the weak-coupling limit there is sensitivity to large 9 signifying a logarithmic divergence, 

U{d)^^^Hd/R), d^a,b. 
For a weakly-coupled sphere in front of a Dirichlet plane, a PFA approximation is 

lOTT d 



(B8) 



(B9) 



(BIO) 
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APPENDIX C: SHORT DISTANCE LIMIT 



1. Cylinders 

In this appendix we want to discuss the short distance limit, for the case of weakly-coupled cylinders, where the 
closest distance between the cylinders is R—a—b = d <^ a,b, which should reduce to the proximity force approximation 
derived in Appendix [B) We will calculate the first correction to the PFA, and compare to the exact result found in 
Sec. IV D] In this limit, we replace the modified Bessel functions by their uniform asymptotic approximants, which in 
leading form yield 

Kl+,n'i^)I^{^a/R)ll,ixb/R) ^ ^ . ttahe-^, (CI) 

on mm [m + m' ) 



where 



and 



i=(l + ^2)-l/2^ i, = (l + z2)-l/2^ i,= (l + z2)-l/2, (C2) 



xa/R xb/R 

Zb = (C3) 



The exponent here is 



m + m' m 



X = 2(m + m')rj{z) — 2mri{za) — 2m''r]{zb), (C4) 



where rj is defined by Eq. (I5.43p . The reason that the force diverges as a + ^ i? is that x vanishes there, for suitable 
values of m and m'. To make this systematic, let us rescale variables, 

m = Ma, m' = M/3, (C5) 

and then when b = R— a, x = when (3a = ab. 

When b = R— a — d, with d small compared to the radius of either cylinder, we assume that the main contribution 
comes from the neighborhood of these values. So we define 

a-^+d, /3=l--^+/3, (C6) 

and we expand the exponent to first order in d and to second order in a and (3 = —a. (The latter constraint ensures 
that a + P = 1.) The result is 

iR a(R — a) 



Then 



which is exactly the result given by the proximity force theorem in Appendix [B] Eq. (jB4p . 

Now we calculate the correction to the PFA. We do this by keeping subleading terms in the uniform asymptotic 
approximation for the product of six Bessel function 



Ki+,n,ix)If^{xa/R)li,ixb/R) - j^-^ 

smnm' m + m' 

uijt 
m + m' 
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where t = t{z) with z = x/(rn + m'), Za = xa/m, zi, = xb/m' , 

3i - 5t^ 



24 



(CIO) 



and X is given by Eq. (jC4[l . Now when we expand x we must go out to order d'*, d?, and d^d. The resuh is 

(PMt 2adMt 



1 



R{R-a) R-a 
dMt{l~t^)R d?Mt^{R-2a)R^ a'^Mt^{l - 2,t^){R^ ~ 'iaR + 2,a^)R^ 



{R-af 3a^iR~a)^ 12a^{R~a)3 

2M2d2rf2^2 ApaH^{R~2aYR^ 2 M^aH^R^ 



{R-a) 



18a4(i? - a) 



3 a^{R-a) 



-{R-2a) 



As above, we replace 



a f R 
m = M— 1 + a — 
it V a 



We expand and ti, in the prefactor using 

dt 

dz 

The PFA is obtained by using the integrals 

dM 



m' = M ( 1 - 



1 -d- 



R 



R-a 



-zf 



dH 

dz"^ 



2t^ - iV'. 



7 = MtR^/a{R-a), 



'2Md/t 



M 



(Cll) 
(C12) 

(C13) 

(C14a) 
(C14b) 



and so, from the expansion we can obtain the result of the integrals over d and M by the algebraic substitutions 

(C15a) 
(C15b) 
(C15c) 



1 

M ~ 


4d t 

-* , Af -> - 

Rt A 


d^ - 


2a{R- a)d 


RH^ ' 


M^d^ - 


a(R — a) 
8Rd ' 



Md 



2-4 



2 1 a(i? — a) 

^ 2 ^ 
3a2(r-a)2 



Md' 



M-'d 



4 3a2(i?-a)2d 
15a3(i?-a)3d 



4 i?4^2 ' 2 

The result is the following correction factor to the PFA in the form given in Eq. (jCSP : 

€ _ R'^ + aR + a'^ d 
U ^ Aa{R-a) R' 



(C16) 



Although this looks slightly different from Eq. (|5.37p . it agrees with the latter when the PFA formula (|5.36p is expressed 
in terms of the form given in Eq. (|C8p . that is, writing d = R — a — b. 



2. Spheres 

Here we see how the proximity force limit is achieved for weakly-coupled spheres. Again, the strategy is to replace the 
modified Bessel functions by their leading uniform asymptotic approximants. The only new element is the a ppe arance 
of the 3-j symbol. Because now only m = appears, there is a very simple approximant for the latter (47l. l48l. l49j : 



/ I' I" 




cos f(Z + + 



2 [{I + + + V - -V + + 1' + 



(C17) 



This result is quite accurate, being within 1% of the true value of the Wigner coefhcient for Vs of order 100 (except 
very near the boundaries of the triangular region, where the approximant diverges weakly). Otherwise, the procedure 
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is rather routine. Letting v = I + 1/2, and similarly for the primed quantities, we expand the exponent resulting from 
the uniform asymptotic expansion about the critical point, with 

v^N{a + a), 1^' ^ N{1- a + a'), i^" ^ N{l + a"), (C18) 

with the constraint a + a' + a" = 0. Replacing the sums over angular momenta by integrals, and changing variables: 

(C19) 



dvdv'dv"^2 dNN^ d(a + d') 



which reflects the restriction emerging from the triangular relation of the Wigner coefficients, a + a' > 0, we find for 
the approximant to the energy when the two spheres are nearly touching: 



AiaA26 2 



4i? TT 7o ^ Jo 



rl-r r°° 1 



a(l — a) 



1/2 



(d + d')i/2 



a \ ^-NtR'(a-a')^ /Aa(R-a) 



\l\2ab 



which is exactly the PFA result (IB9 



Ind, d = R. — a — b <^ a,b, 



(C20) 
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